Nonlinear Mechanical Response of DNA due to Anisotropic Bending Elasticity 
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The response of a short DNA segment to bending is studied, taking into account the anisotropy 
in the bending rigidities caused by the double-helical structure. It is shown that the anisotropy 
introduces an effective nonlinear twist-bend coupling that can lead to the formation of kinks and 
modulations in the curvature and/or in the twist, depending on the values of the elastic constants 
and the imposed deflection angle. The typical wavelength for the modulations, or the distance 
between the neighboring kinks is found to be set by half of the DNA pitch. 
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Genetic information in living cells is carried in the lin- 
ear sequence of nucleotides in DNA. Each molecule of 
DNA contains two complementary strands of nucleotides 
that are matched by the hydrogen bonds of the G-C 
and A-T base pairs (bps). The DNA double- helix can 
be found in several forms that differ from each in the 
geometrical characteristics such as diameter and hand- 
edness. Under normal physiological conditions DNA 
adopts the B-form, in which it consists of two helically 
twisted sugar-phosphate backbones with a diameter 2.4 
nm, which are stuffed with base pairs. The helix is right- 
handed with 10 bps per turn, and the pitch of the helix 
is 3.4 nm. While understanding some aspects of DNA 
functioning requires the details of the chemical bindings 
etc., there are quite a number of processes in cell that 
involve major structural and conformational changes in 
DNA, where typically a small segment of DNA interacts 
with certain proteins |jl|. 

For example, consider the transcription reaction cat- 
alyzed by RNA polymerase, which creates an RNA 
strand that is complementary to the transcribed DNA 
strand. The catalysis takes place inside the transcrip- 
tion bubble: this is an 18 bp stretch of unwound DNA 
helix stabilized by RNA polymerase, which harbors a 12 
bp RNA-DNA hybrid Low resolution structure of E. 
coli RNA polymerase has shown a thumb-like projection 
that seemingly closes around the DNA, and bends it to 
about 90°-120° There are also similar visualizations 
of DNA polymerase, suggesting that the enzyme bends 
a piece of DNA inside it Q. In nucleosomes, a 150 bp 
long segment of DNA wraps about 1.6 turns around the 
histonc complex that has a diameter of about 1 1 nm . 
In all of these DNA-protein complexes, a DNA fragment 
is subjected to a bending of the order of about 4°-7° per 
base pair, which is typically a large curvature. 

The idea of studying the response of DNA to mechan- 
ical stresses is as old as the discovery of the double-helix 
structure itself An elastic description of DNA has 
been developed over the past 20 years, taking on differ- 
ent approaches that include Lagrangian mechanics , 
statistical mechanics M , and molecular dynamics simula- 



tion ||T^. Motivated by the aforementioned structures in 
various DNA-protein complexes, we study the bending of 
short DNA segments using an elastic model. We assume 
that the length-scale set by the helix pitch P ~ 3.4 nm 
that involves 10 base pairs, is long enough so that we can 
use a continuum elastic theory, and yet short enough that 
the spatial anisotropy due to the double-helical structure 
of DNA matters. The elastic model of DNA is used with 
anisotropic bending rigidities, while neglecting the lin- 
ear phenomenological twist-bend ||ll| and twist-stretch 
I p^ couplings. Whereas an isotropic elastic model pre- 
dicts uniform bending and twisting, we find that the 
anisotropy in the bending rigidity introduces an effective 
nonlinear twist-bend coupling [|l^,0 that could lead to 
a variety of bending and twisting structures. 

Depending on the relative strengths of the elastic con- 
stants and the overall amount of deflection, we observe 
six different regimes: (1) A regime in which curvature is 
localized in the form of kinks in a periodic arrangement 
accompanied by simultaneous suppression of twist, while 
the kinks are separated by virtually flat (rod-like) and 
uniformly twisted segments, (2) a regime with synchro- 
nized smooth modulations of the curvature and the twist, 
where a rise in curvature corresponds to suppression of 
twist and vice versa, (3) a regime with uniform curva- 
ture and modulated twist, (4) a regime with the kink-rod 
structure in the curvature and uniform twist, (5) a regime 
with modulated curvature and uniform twist, and (6) a 
regime with both uniform curvature and uniform twist. 
These different "forms" of bent DNA are presented in 
Fig. 1 below, and in Fig. 2 it is specified for what values 
of the parameters each form appears. We find that the 
period of modulations as well as the distance between 
neighboring kinks is always set by half of the DNA pitch, 
i.e. 5 bps, to a good approximation. For a general DNA 
segment bending is always found to reduce the linking 
number due to the effective (nonlinear) twist-bend cou- 
pling. When the length of the DNA segment is equal to 
an integer multiple of 5 bps, however, we observe a tran- 
sition from an over-twisted regime where there are kinks 
located exactly at the end points to an under-twisted 
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regime where the kinks are located exclusively in the in- 
terior. The boundary corresponding to this transition in 
the parameter space of the elastic constants is plotted in 
Fig. 3. 

The elastic model of DNA represents the molecule as 
a slender cylindrical elastic rod. The rod is parametrized 
by the arclength s, and at each point, an orthonormal 
basis is defined with ei, 62, and 63, corresponding to 
the principal axes of the elastic rod. In the undeformed 
state, all the systems are parallel, with the es-axes being 
all parallel to the axis of the rod. The deformation of 
the elastic rod is then parametrized by a mapping that 
relates the local coordinate system to a reference one, 
specified by the Euler angles 0{s), 0(s), and V'(s)- The 
elastic energy is written as 



(3E = 
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ds[Ainl + A2nl + c{n3-ujo?], (i) 



in the dimensionless units, where Ai and A2 are 
the bending rigidities, C is the twist rigidity, and 
ujQ = 2tt/P = 1.85 nm^^ is the spontaneous twist 
of the helix. In terms of Euler angles, we have 
— {d9 /ds) sin ip — (dcj) / ds) sin 6 cos ip , ^l2{s) — 
{d9 / ds) cos i/j+{d(t)/ds) sin 9 sin and 173(3) = {dip/ds) + 
{d(t)/ds) cos 6. Without loss of generality, we assume that 
the rod will stay planar even after the application of a 
bending force, and hence set d(j)/ds — 0. The bending 
is characterized by an overall angle 9q, and we assume 
that the bent segment of DNA is attached to undeformed 
tails, so that it corresponds to bending of a piece of a long 
DNA. Since the tails have a uniform twist characterized 
by Wo, we assume that the twist at the two ends of the 
DNA segment is also llIq, to preserve continuity. 

After defining A = ^{Ai + A2), A' = \{A2 ~ Ai), and 
setting s = Lt, we have 



f3E 



2L 



dt 



C 



A'9l 



(2) 



The above 



where we have used 9 = ^ and ij} = 
equation shows manifestly that there are three effective 
dimensionless parameters in the system: A/A' , C/{A'9'^), 
and loqL. We can find the extremum of [3E by applying 
standard variational techniques to Eq. (^. The resulting 
Euler-Lagrange equations for 9{t) and i}{t) are found as 



9 _ 2ijj{9/9o)sin2'iP 
9^> " (^+cos2^) ' 



(3) 



The above equations should be solved subject to the 
boundary conditions '0(0) = ip(l) = (^qL, 9{0) = and 
9{1) — 9q. While the mean bending rigidity A is known 
to be about 50 nm in the salt-saturation limit 1^, and 
the twist rigidity C is believed to be roughly 75 nm , 
little is known about the parameter A' that describes 
the anisotropy. Hence, we solve the above set of cou- 
pled nonlinear equations numerically for various values 
of the parameters A /A' and C/{A'9'^), and for several 
DNA lengths. 
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FIG. 1. Different possible conformations of the DNA for 
various values of C/(A'e»o) and A/ A'. The DNA length is 
18 bps. The solid lines show the normalized curvature 9/6o 
and the dotted lines show the twist 7/). (a) Kink-rod cur- 
vature structure and a corresponding singular suppression of 
twist at the kinks. The plot corresponds to A/ A' — 1.02 and 
C/{A'eo) = 0.5. (b) Modulated curvature and twist. The 
plot corresponds to A/A' = 3 and C/iA'Sg) = 0.5. (c) Uni- 
form curvature and modulated twist. The plot corresponds to 
A/ A' = 79 and C/{A'9o) = 0.5. (d) Kink-rod curvature struc- 
ture and uniform twist. The plot corresponds to A/A' = 1.02 
and C/{A'6q) — 30. (e) Modulated curvature and uniform 
twist. The plot corresponds to A/A' = 3 and C/(A'6'o) = 30. 
(f) Uniform curvature and uniform twist. The plot corre- 
sponds to A/A' = 79 and C/{A'e'^) = 30. 
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(4) 



We observe that the above equations can acquire more 
than one solution depending on the values of the parame- 
ters, with the different DNA conformations characterized 
by deferent energies and different values of the linking 
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number, defined as Lk = [tp{l) - ?/;(0)]/(27r) jl^. Note 
that for an undeformed DNA we have Lkg = ujoL/{2tt). 
For each value of DNA length that is not an integer mul- 
tiple of 5 bps, we find six distinct types of behaviors 
for various values of A/ A' and C/{A'9q), as presented in 
Fig. 1 and described above. For all of these conforma- 
tions, which correspond to the lowest energy case among 
the different possible solutions, the change in the linking 
number ALk = Lk — Lko is negative, and DNA is thus 
under-twisted. 

We also observe in all of these cases that the peri- 
odicity of the curves, whether there are modulations or 
kinks, is set by half of the DNA pitch, i.e. 5 bps, to a 
good approximation. This can be understood by looking 
at the first-integrals of Eqs. (||) and which can be 
calculated as 



9o (-^+cos2V') 



and 



(:4%) 



COS 2-0) 



a2, 



(5) 



(6) 



respectively, where ai and a2 are the integration con- 
stants. Using a linear approximation ■0 ~ 27rs /P, we can 
then expect a periodic behavior for both 9 and in Eqs. 
(|) and d), with a period of P/2. 
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FIG. 2. A diagram delineating the different regimes corre- 
sponding to the six distinct forms of bent DNA, depending 
on the values of C/{A'e^) and A/ A'. The DNA length is 
18 bps. For fixed values of the elastic constants A, A' , and 
C, bending of DNA at various angles would correspond to 
vertical lines. Note that in this case, the only permissible 
transitions are from d-form to a- form, from e-form to b-form, 
and from f-form to c-form, because the zone boundaries are 
also vertical. 



When the ratio A/ A' (that is by definition always 
greater than 1) is somewhat close to 1, the combination 
A/ A' + cos 2?/' could be close to zero, and then Eq. (||) 
implies that 9/9q must be large. This corresponds to 
the kinks in the curvature. In this case, for sufficiently 
small C/{A'9l) (numerically we find < 1.6), Eq. (||) im- 
plies that the same effect should simultaneously happen 
in the twist, which corresponds to Fig. la. For larger 
C/{A'9'^), however, this pattern should disappear from 
the twist structure leaving only a uniform twist, which 
corresponds to Fig. Id. 

For intermediate values of A/A\ the periodicity im- 
posed by the term A/ A' -f cos 2il> is still felt in the curva- 
ture, in the form of modulations. Again, in this case, a 
small value for C/{A'9l) leads to accompanying modula- 
tions in the twist, corresponding to Fig. lb, while large 
C / {A'9q) will wash it out to uniform twist, corresponding 
to Fig. le. Note that an increase in the curvature always 
coincides with a decrease in the twist, and vice versa, as 
can be manifestly understood from Eqs. (||) and (^). 

When A/ A' is very large, the dominant term in the 
denominator of Eq. ^ is A/A', and hence the curva- 
ture approaches a constant value. Since the parameters 
are normalized such that a constant curvature should be 
equal to 1, we find that in this limit ai ~ A/A'. Then 
a small value of C/{A'9f^) can make it possible for the 
modulations to be felt in the twist, corresponding to Fig. 
Ic, whereas for large values of C/{A'9q) these modula- 
tions are also washed out and we have a uniform twist, 
which corresponds to Fig. If. 

In Fig. 2, we have summarized the above results in a 
diagram that delineates the different regimes correspond- 
ing to the different forms of bent DNA defined in Fig. 1, 
for various values of C /{A'Oq) and A/A'. Note that while 
the horizontal line (the boundary that separates a, b, and 
c regions from d, e, and f regions) is slightly tilted, the 
other two lines are virtually vertical within our numerical 
precision. Since for fixed values of the elastic constants 
A, A', and C, bending of DNA at various angles would 
correspond to vertical lines, the fact that the two zone 
boundaries are vertical requires that the only permissi- 
ble transitions upon bending are from d-form to a-form, 
from e-form to b-form, and from f-form to c-form. This 
could provide an interesting possible experimental check 
of the present analysis. 

In the a-form regime, when the DNA length is between 
5?! bps and 5{n+ 1) bps with n being an integer number, 
the number of kinks is rt + 1. One could then ask what 
happens if the DNA length is exactly equal to 5n bps. 
In this case, two different sets of solutions are observed, 
as shown in the insets of Fig. 3: (1) A set of solutions 
in which there are n kinks situated exclusively in the in- 
terior, leaving an almost vanishing curvature at the end- 
points. The value for ALk in this case is negative, which 
means that DNA is under- twisted. (2) Another set of 
solutions in which there are n — 1 kinks in the interior, 
and the other kink is splitted in half and located exactly 
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at the end-points, leading to a large curvature at the two 
ends. In this case, ALk is positive, i.e. DNA is over- 
twisted. We found that both of these sets of solutions 
could correspond to the minimum energy case, depend- 
ing on the values of the parameters. The transition line 
between the two regimes is plotted in Fig. 3, for a DNA 
segment of 15 bps. 

In the present treatment thermal fluctuations have 
been neglected. This is justified by noting that the typ- 
ical lengths considered here are much smaller than the 
persistence length of DNA, which is about 150 bps, and 
therefore we do not expect considerable deformation fluc- 
tuations of the DNA. However, due to the presence of 
multiple energy minima that are not too different from 
each other, we do expect that the conformation of DNA 
jumps between the different possible minimum energy so- 
lutions, provided that the barrier is not too large. This 
could be particularly interesting if the bent DNA is in- 
volved in a dynamical process, such as transcription [n8|. 



L = 15 bps 




FIG. 3. A diagram showing the different conformations of 
DNA, when its length is an integer multiple of 5 bps. A tran- 
sition is observed from an over-twisted regime where there are 
kinks located exactly at the end-points, to an under-twisted 
regime where the kinks are located exclusively in the interior. 

In conclusion, we have studied the response of a fi- 
nite DNA segment to bending, taking into account the 
anisotropy in the bending rigidity caused by the double- 
helical structure. We have shown that the nonlinear 
twist-bend coupling caused by the anisotropy can lead 
to formation of kink-rod patterns and modulations in the 
curvature and/or twist. This effect may be related to the 
recent observation in high resolution X-ray, that DNA is 
not uniformly bent around the histone octamer and its 
curvature is larger at some places in the nucleosome 
We finally note that the present analysis can be also ap- 
plied to other biopolymers — the best candidate for an 



experimental realization of these effects will probably be 
F-actin, which has a much longer persistence length and 
can thus be more easily manipulated in its stiff limit. 

We are grateful to M.A. Jalali and T.B. Liverpool for 
helpful discussions and comments. 
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